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We study the scattering of light pseudoscalar mesons (tt, K) off charmed mesons (D, 
D s ) in full lattice QCD. The S-wave scattering lengths are calculated using Liischer's finite 
volume technique. We use a relativistic formulation for the charm quark. For the light 
quark, we use domain-wall fcrmions in the valence sector and improved Kogut-Susskind sea 
quarks. We calculate the scattering lengths of isospin-3/2 Dir, D s tt, D s K, isospin-0 DK and 
isospin-1 DK channels on the lattice. For the chiral extrapolation, we use a chiral unitary 
approach to next-to-lcading order, which at the same time allows us to give predictions for 
other channels. It turns out that the interpretation of the £>* (2317) as a DK molecule is 
consistent with the results. At the same time, we also update a prediction for the isospin 
breaking hadronic decay width r(D* (2317) D s ir) to (89 ± 27) keV. 
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I. INTRODUCTION 

In 2003, the BaBar Collaboration discovered a positive-parity scalar charm-strange meson 
Z)* (2317) with a very narrow width [1]. The state was confirmed later by the CLEO Collab- 
oration [2]. The discovery of this state has inspired heated discussions in the past decade. The 
key point is to understand the low mass of this state, which is more than 100 MeV lower than 
the prediction for the lowest scalar cs meson in, for instance, the Godfrey- Isgur quark model [3]. 
There are several interpretations of its structure, such as being a DK molecule, the chiral partner 
of the pseudoscalar D s , a conventional cs state, coupled-channel effects between the cs state and 
DK continuum etc. For a detailed review of the properties and the phenomenology of these states 
see Ref. [4] . In order to distinguish them, one has to explore the consequences of each interpre- 
tation, and identify quantities which have different values in different interpretations. Arguably 
the most promising quantity is the isospin breaking width r(D* (2317) — > D s tt). It is of order 
10 keV if the £>* (2317) is a cs meson [5, 6], while it is of order 100 keV [7-9] in the DK molecular 
picture due to its large coupling to DK and the proximity of the DK threshold. Thus, the study 
of DK interaction is very important in order to understand the structure of D* (2317). Although 
a direct simulation of the DK(I = 0) channel suffers from disconnected diagrams, one may obtain 
useful information on the DK interaction by calculating the scattering lengths of the disconnected- 
diagram-free channels which can be related to DK{1 = 0) through SU(3) flavor symmetry. This 
is the strategy we will follow in this paper. 

Lattice QCD calculations of the properties of hadronic interactions such as elastic scattering 
phases shifts and scattering lengths have recently started to develop. Precision results have been 
obtained in the light meson sector for certain processes such as pion-pion, kaon-kaon and pion- 
kaon scattering and preliminary results for baryon-baryon scattering lengths have been presented. 
A review of these calculations can be found in Ref. [10]. In the heavy meson sector, only a few 
calculations have been done, including quenched calculations in Refs. [11, 12] and calculations in 
full QCD in Refs. [13, 14]. In this work, we study scattering processes where one of the hadrons 
contains a charm quark in full lattice QCD. 

Extracting hadronic interactions from Lattice QCD calculations is not straightforward due to 
the Maiani- Testa theorem [15], which states that the S'-matrix can not be extracted from infinite- 
volume Euclidean-space Green functions except at kinematic thresholds. However, this problem 
can be evaded by computing the correlation functions in a finite volume. Liischer has shown that 
one can obtain the scattering amplitude from the energy of two particles in a finite volume [16, 17]. 
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We use Liischer's finite volume technique to calculate the scattering lengths. We then use unitarized 
chiral perturbation theory to extrapolate our results to the physical pion mass. Having fitted the 
appearing low-energy constants (LECs) to the lattice data, we are also able to make predictions for 
other channels, in particular for the isospin zero, strangeness one channel in which the Z)* (2317) 
resides. 

The paper is organized as follows. The lattice formulation of the light and heavy quark ac- 
tions will be discussed in Section II. Liischer's formula will be briefly introduced in Section III. 
The numerical results for the scattering lengths of five channels DK(I = 0), DK(I = 1), D S K, 
Dir(I = 3/2) and D s tt which are free of disconnected diagrams will be given in Section IV. Chiral 
extrapolations will be performed using unitarized chiral perturbation theory in Section V, and val- 
ues of the LECs in the chiral Lagrangian will be determined. Predictions for other channels using 
the LECs are given Section VI, and in particular, implications on the L>* (2317) will be discussed. 
The last section is devoted to a brief summary. 

II. LATTICE FORMULATION 
A. Light-Quark Action 

In this work we employ the "coarse" (a ~ 0.125 fm) gauge configurations generated by the 
MILC Collaboration [18] using the one-loop tadpole-improved gauge action [19], where both 0(a 2 ) 
and 0(g 2 a 2 ) errors are removed. For the fermions in the vacuum, the asqtad-improved Kogut- 
Susskind (staggered) action [20-25] is used. This is the so-called Naik action [26] {0(a 2 ) improved 
Kogut-Susskind action) with smeared links for the one-link terms so that couplings to gluons with 
any of their momentum components equal to ir/a are set to zero, resulting in a reduction of the 
flavor symmetry violations present in the Kogut-Susskind action. 

For the valence light quarks (up, down and strange) we use the five-dimensional Shamir [27, 28] 
domain- wall fermion action [29]. The domain- wall fermion action introduces a fifth dimension of 
extent L5 and a mass parameter M5; in our case, the values L5 = 16 and M5 = 1.7, both in lattice 
units, were chosen. The physical quark fields, q(x,t), reside on the 4-dimensional boundaries of 
the fifth coordinate. The left and right chiral components are separated on the corresponding 
boundaries, resulting in an action with chiral symmetry at finite lattice spacing as L5 — > 00. We 
use hypercubic-smeared gauge links [30-33] to minimize the residual chiral symmetry breaking, 
and the bare quark-mass parameter (am)J wf is introduced as a direct coupling of the boundary 
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Ensemble 


P 


ami 


am s 


amf wt 




A'cfgs 


1 v props 


m007 


6.76 


0.007 


0.050 


0.0081 


0.081 


461 


2766 


mOlO 


6.76 


0.010 


0.050 


0.0138 


0.081 


636 


3816 


m020 


6.79 


0.020 


0.050 


0.0313 


0.081 


480 


1920 


m030 


6.81 


0.030 


0.050 


0.0478 


0.081 


563 


1689 



TABLE I: The parameters of the configurations and domain-wall propagators used in this work. The 
subscript I denotes light quark, and s denotes the strange quark. The superscript "dwf" denotes domain- 
wall fermion. 

chiral components. The light quark propagators were provided to us by the NPLQCD [10] and 
LHP [34-36] Collaborations. 

The calculation we have performed, because the valence and sea quark actions are different, 
is inherently partially quenched and therefore violates unitarity. Unlike conventional partially 
quenched calculations, to restore unitarity, one must take the continuum limit in addition to 
tuning the valence and sea quark masses to be degenerate. This process is aided by the use of 
mixed-action chiral perturbation theory [37-42]. Given the situation, there is an ambiguity in the 
choice of the valence light-quark masses. One appealing choice is to tune the valence light quark 
masses such that the valence pion mass is degenerate with the Goldstone staggered pion mass. 
In the continuum limit, the Nf = 2 staggered action has an SU(8)l <8> SU(8)r U(l)y chiral 
symmetry due to the four-fold taste degeneracy of each flavor, and each pion has 15 degenerate 
partners. At finite lattice spacing this symmetry is broken and the taste multiplets are no longer 
degenerate, but have splittings that are O(a 2 s o?) [20-22, 25, 43]. The propagators used in this work 
were tuned to give valence pions that match the Goldstone Kogut-Susskind pion. This is the only 
pion that becomes massless in the chiral limit at finite lattice spacing. As a result of this choice, 
the valence pions are as light as possible, while being tuned to one of the staggered pion masses, 
providing better convergence in the chiral perturbation theory needed to extrapolate the lattice 
results to the physical quark-mass point. This set of parameters, listed in Table I, was first used 
by LHPC [34, 35] and recently utilized to compute the spectroscopy of hadrons composed of up, 
down and strange quarks [36]. 
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B. Heavy-Quark Action 

For the charm quark we use a relativistic heavy quark action motivated by the Fermilab ap- 
proach [44]. This action controls discretization errors of 0{{am,Q) n ). Following the Symanzik 
improvement [45], an effective continuum action is constructed using operators that are invariant 
under discrete rotations, parity-reversal and charge-conjugation transformations, representing the 
long-distance limit of our lattice theory, including leading finite-a errors. Using only the Dirac 
operator and the gluon field tensor (and distinguishing between the time and space components of 
each), we enumerate seven operators with dimension up to five. By applying the isospectral trans- 
formations [46] , the redundant operators are identified and their coefficients are set to appropriate 
convenient values. The lattice action then takes the form 



S = So + S B + S E , (1) 



with 



S = ^2Q(i 



to 



+ (70 V - ^A ) + (j* V * " | Ai ) 



Q(x), (2) 



i<j 



Se = -^c e ^ Q(x) ^E a oiF }j Q(x) , 



(4) 



where the operator Q(x) annihilates a heavy quark field, a is the lattice spacing, Vo and V, are 
first-order lattice derivatives in the time and space directions, Ao and Aj are second-order lattice 
derivatives, and F^ v is the gauge field strength tensor. The spectrum of heavy-quark bound states 
can be determined accurately through \p\a and (otoq)™ for arbitrary exponent n by using a lattice 
action containing too, v, cb and ce, which are functions of amg. 

The coefficients cb and ce are different due to the broken space-time interchange symme- 
try, which can be computed in perturbation theory by requiring elimination of the heavy-quark 
discretization errors at a given order in the strong coupling constant a s . We use the tree- level 
tadpole-improved results obtained by using field transformation (as in Ref. [46]): 

cb = 4, c E = hl + v)^, (5) 

where no is the tadpole factor 

1/4 



U0= \3 

\ v 



W^{U P )\ , (6) 
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and Up is the product of gauge links around the fundamental lattice plaquette p. The remaining 
two parameters mo and u are determined nonperturbatively. The bare charm-quark mass mo is 
tuned so that the experimentally observed spin average of the J/tp and r\ c masses 



is reproduced. For each ensemble, we calculate M avg at two charm-quark masses (denoted 
m\ = 0.2034 and 777.2 = 0.2100) and linearly extrapolate it to the experimental value to deter- 
mine the parameter mo = rric hys . The value of v must be tuned to restore the dispersion relation 
E\ = m\ + c 2 p 2 such that c 2 = 1. To do this, we calculate the single-particle energy of rj c , 
J/ip, D s and D at the six lowest momenta (with unit of a -1 ): (2ir/L)(0, 0, 0), (27r/L)(l, 0, 0), 
(27r/L)(l,l,0),(27r/L)(l,l,l),(27r/L)(2,0,0),(2vr/L)(2,l,0). For each ensemble, the energy lev- 
els are calculated at the two charm-quark masses (mi and 771,2 ) and extrapolated to the physical 
charm-quark mass mc hys . The values of c 2 are obtained by fitting the extrapolated energy levels 
to the dispersion relation. We tune v using the dispersion relation of r\ c . The dispersion relations 
for either the charmonium J/ip or the charm- light mesons (D and D s ) are generally consistent 
with c 2 = 1 within 1-2%. Since the values of v and mo are coupled, one needs to iterate the 
tuning process in order to achieve a consistent pair of values. For the details of tuning the bare 
charm-quark mass uiq and the value of v, see reference [47]. 



Liischer has shown that the scattering phase shift is related to the energy shift (AE) in the 
total energy of two interacting hadrons in a finite box [16, 17]. 

The center-of-mass momentum p can be obtained by the relation 




(7) 



III. 



LUSCHER'S FORMULA 




(8) 



where mi and m2 are the rest masses of the two hadrons. 



To obtain pcot5(p), where S(p) is the phase shift, we use the formula [48] 




(9) 



where the S function is defined as 




UK a 



(10) 
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The sum is over all three- vectors of integers j such that |j| < A, and the limit A — > oo is implicit. 
If the interaction range is much smaller than the lattice size, p cot S(p) is given by 

pcotS(p) = ^ + 0(p 2 ), (11) 

where a is the scattering length (not to be confused with the lattice spacing which has the same 
notation and dimension). Note that we take the sign convention that a repulsive interaction has a 
negative scattering length. The higher order terms in Eq. (11) can be ignored if the effective range 
of the interaction is much smaller than the length scale associated to the center-of-mass momentum 
p. If we ignore the higher order terms, the scattering length can be calculated by 

— "-'((gn- < i2 » 

IV. NUMERICAL RESULTS 

In the following, we list all the channels we study. The interpolating operators for these two 
particle states are 

O d , v = DJk + , I = t 3/2 = D + tt + , DsK = DfK+, 
°dR = D + K°, O 1 ^ = D+K- - D°K°, 

where D~, Df, D + , K°, K~, K + and ir + are the operators for one particle states, the subscripts 
7r, D, K and K represent the isospin triplet (ir + , ir°, ir~) and doublets (D + , D°), (K + , K°) and 
(K°, K~), respectively. 

The total energy of two interacting hadrons (h\ and /12) is obtained from the four-point corre- 
lation function: 

G h ^(t) = (O hlh2 (tyO hlh2 (0))- (13) 
To extract the energy shift AE, we define a ratio R hlh2 (t): 

^M- cC^M ^"*-^ (14) 

where G hl (t,0) and G h2 (t, 0) are two-point functions. AE is obtained by fitting R hlh2 (t) to a 
single exponential in the region where the effective mass exhibits a plateau. 

For each channel, we calculate the ratio R hlh ' 2 at two different charm quark masses and four 
different light valence quark masses. Figure 1 shows the effective energy shifts of each channel 
calculated from ensemble m007 at the bare charm-quark mass 1712 = 0.2100. The fitted energy 
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shifts and the fitting ranges are indicated by the grey bars in these plots. The heights of the grey 
bars show the statistical errors. The effective energy shift plots for other ensembles are similar. 
The energy shifts are linearly extrapolated to the physical charm-quark mass on each ensemble. 
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FIG. 1: Effective energy shifts plots of the scattering channels D s ir, D S K, DK(I = 0), DK(I = 1), 
Dtt(I = 3/2). All plots are for ensemble m007. The grey bars show the fitted energy shifts and the fitting 
ranges. The height of the grey bars show the statistical errors. 
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V. CHIRAL EXTRAPOLATIONS OF THE SCATTERING LENGTHS 

Because the simulations are performed at unphysical quark masses, chiral extrapolation is nec- 
essary in order to obtain the values of scattering lengths at the physical quark masses. There 
have been calculations based on chiral Lagrangians for these scattering lengths [49-52] . They were 
first calculated in Ref. [49] using a unitarized chiral approach. The basic observation is because 
of the coupled-channel effect and the large kaon mass, the interaction of some of the channels is 
so strong that a nonperturbative treatment is necessary, and in one channel even a bound state 
is produced. The method was followed up recently in Ref. [52]. Some other authors treated the 
interaction perturbatively, and calculated the scattering lengths up to leading one-loop order in 
chiral perturbation theory with [50] and without [51] a heavy quark expansion. Here we take the 
same route as Ref. [49], and resum the chiral amplitude up to next-to- leading order, which is 0(p 2 ). 
The resummed amplitude in the on-shell approximation reads [53-55] 

T(s) = V(s)[l-G(s)V(s)}-\ (15) 

where V(s) is the S'-wave projection of the 0(p 2 ) scattering amplitude, and G(s) is the scalar loop 

function regularized by a subtraction constant a(A) 

m s 1 f~/,x , m \ m 2 — m\ + s^m 2 c r 2 2 , 

= 16^r (A) + ln + 2s ln m\ + 2s N* + + 

— ln(— s + m\ — m\ + a) + ln(s + m\ — m\ + a) — ln(— s — m\ + m\ + a)] |, (16) 

with a = {[s — (mi + m,2) 2 ][s — (mi — m^) 2 ]} 1 ^ 2 . A is the scale of dimensional regularization, and 
a change of A can be absorbed by a corresponding change of a(A). The value A = 1 GeV will be 
taken in the following. Promoting T(s), V(s) and G(s) to be matrix- valued quantities, it is easy 
to generalize Eq. (15) to coupled channels. 

Using the 0(p 2 ) chiral Lagrangian constructed in Ref. [9], the scattering amplitudes are given 

by 

1 



V(s,t,u) - p2 



' LO (a - u) - 4Q)/io + 2Cihi - 2C 2i H 24 (s, t, u) + 2C 35 #35 (s, t, u) 



(17) 



4 

where F is the pion decay constant in the chiral limit, and the coefficients Cj can be found in 
Table II. Further, 

H 2 i(s, t, u) = 2h 2 p2 ■ Pa + h A (pi ■ p 2 P3 -Pi+Pi- PiP2 ■ P3), 

and 

H 35 (s, t, u) = h 3 p 2 ■ p 4 + h 5 (pi ■ p 2 p 3 -pa+Pi- PiP2 ■ Pz)- 
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TABLE II: The coefficients in the scattering amplitudes V(s, t, u). The channels arc labelled by strangeness 
(5) and isospin (7). 

Note that the term h\\+ = h\(\+ — (x+)/3) in the Lagrangian in Refs. [9, 49] has been replaced 
by hix+, which amounts to a redefinition of ho (for the details of the Lagrangian and the definition 
of x+, we refer to Refs. [9, 49]). In this way, the h\ term does not contain the 1/N C , with N c being 
the number of colors, suppressed part (%+) any more. This was also done in Ref. [52]. 

In previous works [9, 49], the large- N c suppressed low-energy constants (LECs) ^0,2,4 were 
dropped to reduce the number of parameters. However, when fitting to the lattice data at several 
unphysical quark masses, this is no longer necessary. In this work, we will keep all of the LECs 
at this order, and fit them to the lattice data. This strategy were also taken in Refs. [50-52], 
where the preliminary lattice results [13] were used. By definition, the LECs are independent of 
the pion mass. We further need to assume that the subtraction constant is the same for various 
channels, and neglect its pion mass dependence. In principle, this assumption is not necessary for 
a unitarization procedure matched to the full one-loop level of the perturbative calculation [55, 56], 
which will be left for the future. 

From the SU(3) mass splitting of the charmed mesons, the value of h\ is fixed to be h\ = 0.42. 
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DK(I = 1) 


DK(I = 0) 


D S K 


L>tt(/ = 3/2) 


D s ir 


M n 


a (fm) 


m007 


-1.19(0.40) 


5.34(1.45) 


-1.58(0.14) 


-1.16(0.30) 


0.08(0.04) 


0.1842 


0.1207 


mOlO 


-1.89(0.12) 


6.21(1.04) 


-1.55(0.09) 


-1.38(0.10) 


0.08(0.03) 


0.2238 


0.1214 


m020 


-1.49(0.25) 


4.43(1.33) 


-1.40(0.20) 


-1.08(0.30) 


0.13(0.05) 


0.3113 


0.1202 


m030 


-1.59(0.13) 


7.46(1.56) 


-1.67(0.10) 


-1.68(0.13) 


0.32(0.05) 


0.3752 


0.1200 



TABLE III: The values of scattering lengths for five channels and in lattice units. The values of the 
lattice spacing a are also given in the last column. 

We still have six parameters, which are a, /13, /15, ho, hi and /14. They are to be fitted to the lattice 
data. However, there is a high correlation between /13 and /is, as well as a similar correlation 
between /12 and /14. In the heavy quark limit, the S*-wave projected amplitudes cannot distinguish 
the /i4( 5 ) terms from the /i 2 (3) ones [57]. Hence, we may reduce the correlations largely by rewriting 
H24,(s,t,u) and H^(s,t,u) as 

H 2i (s, t, u) = 2h 2 4P2 -PA + h A (pi ■ p 2 P3 -Pi+Pi- PaP2 ■ P3 ~ 2M|,P2 • Pa) , 

and 

H 3 b{s, t, u) = h 35 p 2 -p4 + h 5 (pi ■ p 2 P3 -P4+Pi- PaP2 ■ P3 ~ 2Mpp 2 ■ pa) , 

where M D = (M^ + M]J y )/2, the average of the masses of the D and D s at the physical pion 
mass, is introduced to match the dimensions. The new parameters h,2A and /135 are dimensionless, 
and their relations to the old ones are /124 = /12 + h' A and /135 = h% + 2h' 5 , where h' 4 = h^M^ and 
h' 5 = h b M 2 D . 

There are four different light quark masses in our data set, corresponding to the four ensembles 
(m007, mOlO, m020 and m030) with pion masses approximately 301 MeV, 364 MeV, 511 MeV and 
617 MeV, respectively There are in total 20 data points in the five channels. The values of the 
scattering lengths for all the channels are collected in Table III, where the values of M n and lattice 
spacings for all the ensembles are also given. However, for a pion mass as high as 617 MeV, the 
kaon mass would be even higher, around 700 MeV. Such values are too large for a controlled 
chiral expansion. Therefore, we will only fit to the ensembles m007, mOlO and m020. To minimize 
the contamination from a particular scale-setting method, we fit to the dimensionless product of 
the pion mass and the scattering length. The fit was performed using the FORTRAN package 
MINUIT [58]. The best fit has x 2 /dof = 1.17, and the resulting parameters are collected in 
Table IV, where the asymmetric la uncertainties are calculated using the MINOS algorithm in 
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TABLE IV: Results of fitting to the lattice data of the scattering lengths with 6 parameters. 
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FIG. 2: Fit to the data of the scattering lengths corresponding to ensembles m007-m020 in each channel. 
The superscript (5, /) is the (strangeness, isospin) for each channel. 



MINUIT. A comparison of the best fit and the lattice data is shown in Fig. 2, where the solid 
curves correspond to the results of the best fit, and the bands reflect the uncertainties propagated 
from the lattice data. At the physical pion mass, the extrapolated scattering lengths for the five 
channels are presented in Table V. 

One can see that all the dimensionless parameters have a natural size, i.e., the absolute values 
of /io,24,35 and h' 3 5 are of order unity. However, with the current fit, some of the parameters are 
not well constrained. For instance, one of the asymmetric uncertainties of either ho or /124 is 
much larger than its central value. It turns out that when ho takes a larger value, /i24 will take a 
smaller value, and a will also take a smaller value. The asymmetry in the uncertainties of these 
three parameters is larger than that for the others. This indicates a strong correlation of these 
parameters. Furthermore, for a large part of the uncertainty range, h,Q is larger than h\. This is in 
contrast to the N c counting which suggests an opposite ordering. As will be shown in Section VI B, 
the parameters can be much better constrained with an additional assumption. 

When performing the fit, we have used the physical value for the pion decay constant F = 
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Channels 


DK(I = 1) DK(I = 0) D S K Dir(I = 3/2) D s ir 


a (fm) 


_ 21 +0.02 oq0 +0.20 _ 02f) +0.02 _ n 10 o+0.005 -0 002+ 010 
u - zl -0.01 u - yu -0.32 u ' zu -0.01 u - iuo -0.002 U ' UUZ -0.002 


TABLE V: The scattering lengths extrapolated to the physical light quark masses. 


Channels 


Dn(I = 1/2) DK(I = 1) D S K 


a (fm) 


0.3llg;g 0.09™ + i(0.15lg : «) -0.06±S:8g + <(0.47±g:S) 



TABLE VI: Scattering lengths of Dn(I = 1/2), DK(I = 1) and D S K at the physical pion mass predicted 
from the fit. 

92.21 MeV [59]. The difference from the chiral limit value and hence its pion mass dependence is 
a higher order effect, and is neglected here, although it might have some influence. Furthermore, 
we have expressed the masses of the other hadrons in terms of masses in the SU(2) chiral limit, 

° 9 ° 

and the pion mass. The kaon mass is written as Mk = Mk + M* / (AM k) (see, for instance, [60]), 
and the charmed meson masses are 

M 2 ° M 2 

M D = M D + (hi + 2h )^, M Ds = M Ds + 2h . (18) 

M D M Ds 

VI. IMPLICATIONS FOR OTHER CHANNELS 
A. Scattering lengths 

In this work we did not calculate the scattering lengths on the lattice of the channels whose 
Wick contractions involve disconnected diagrams due to the computational difficulties, as well as the 
additional lattice artifacts present in these channels due to the use of Kogut-Susskind sea quarks. 
However, once we have determined the LECs in the chiral Lagrangian, we can make predictions 
on the scattering lengths of these channels. The results for the scattering lengths of Dtt(I = 1/2), 
DK(I = 1) and D S K at the physical pion mass are presented in Table VI. For DK(I = 1), the 
imaginary part of the scattering length originates because it couples to D s ir with a lower threshold. 
Similarly, D S K couples to Dtt and Drj so that the scattering length is complex, too. The result for 
the Dtt(I = 1/2) channel is consistent with the indirect extraction from lattice calculations of the 
Dtt scalar form factor (0.41 ± 0.06) fm [61]. At a pion mass of about 266 MeV, our prediction is 
(0.90+°;^) fm, compatible with the very recent full QCD calculation (0.81 ± 0.14) fm [14]. 
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Fitting range 


X 2 /dof 


h 


h24 


K 


^35 




m007-m020 


1.20 


-0 09+ - 08 


-o ofi+ 003 

u - uo -0.04 


_1 cq+o.60 
i.oy_ 77 


U ' 10 -0.13 


-3.n±i:SS 



TABLE VII: Results of fitting to the lattice data of the scattering lengths with 5 parameters. The subtraction 
constant is solved from fixing the pole in the (S, I) = (1, 0) channel to 2317.8 MeV. 

The most interesting channel is the one with (S,I) = (1,0), where the D* (2317) resides. 
This state was proposed to be a hadronic molecule with a dominant DK component by several 
groups [62-66]. The attraction in this channel is so strong that a pole emerges in the resummed 
amplitude. However, as can be seen from Table IV, the uncertainties of some of the parameters 
are rather large, which causes a large range of the pole position. The pole can be a bound state 
pole in the first Ricmann sheet at around 2200 MeV, and it can also be a resonance pole in the 
second Riemann sheet. As a result, the scattering length in this channel ranges from — oo to +oo. 
A large value of positive scattering length implies the existence of a virtual state, which is a pole 
on the real axis on the second Riemann sheet, while a large value of negative scattering length (the 
absolute value is large) implies a bound state close to threshold. As emphasized in, for instance, 
Refs. [67, 68], if there is an S- wave shallow bound state, the scattering length is related to the 
binding energy, and to the wave function renormalization constant Z, with (1 — Z) being the 
probability of finding the molecular component in the physical state (for Z = 0, the physical state 
is purely a bound state). The relation reads 




where \i and e are the reduced mass and binding energy, respectively. Corrections of the above 
equation come from neglecting the range of forces, 1//3, which contains information of the D s r\ 
channel. Were the D* (2317) a pure DK bound state, the value of DK(I = 0) scattering length 
would be a = —1.05 fm. 

B. Isospin breaking width of the D* (2317) 

In the following, we will assume that the D* (2317) corresponds to the pole generated in the 
(S,I) = (1,0) channel, and explore the implications of our lattice calculation on this state. We 
will fix the pole position to the mass of the D* (2317), 2317.8 MeV [59], on the first Riemann 
sheet. We fit the lattice results of the scattering lengths with ho, /124, ^35, h' 4 and h' 5 , and adjust 
the subtraction constant a(X = 1 GeV) to reproduce the mass of the D* (2317). Again, we only 
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FIG. 3: Fit to the data corresponding to ensembles m007-m020 in each channel with 5 parameters. The 
subtraction constant is solved from fixing the pole in the (S,I) = (1,0) channel to 2317.8 MeV. 



Channels 


DK (I = 1) 


DK (I = 0) 


D S K 


Dn{I = 3/2) 


D S TT 


a(fm) 


-0.20(1) 


0.77(12) 


-0.18(1) 


-0.099(2) 


-0.002(1) 



TABLE VIII: The scattering lengths extrapolated to the physical light quark masses from the 5-parameter 
fit. 

fit to the ensembles m007, mOlO and m020. The best fit gives x 2 /dof = 1.20, which is almost the 
same as the one with one more parameter in Section V. The parameter values together with the 
la statistical uncertainties are given in Table VII. The fitted results are presented in Fig. 3. At 
the physical pion mass, the extrapolated values of the scattering lengths are listed in Table VIII. 

With the newly fitted parameters, the scattering lengths for several other channels are predicted, 
and the results are listed in Table IX. One sees that the value for the DK{I = 0) channel is close 
to the result of Eq. (19), —1.05 fm, with Z = 0. The derivation from this value is partly due to the 



Channels 


Dn(I = 1/2) 


DK(I 


= 0) 


DK(I = 1) 


D S K 


a (fm) 


0.37 ±0.01 


-0.85: 


-0.07 
-0.05 


°- 03 -o:o5 +«(0.05±0.04) 


-0.19 ± 0.05 + «(0.12±g;g|) 



TABLE IX: Scattering lengths of Dtt(I = 1/2), DK(I = 0), DK(I = 1) and D S K at the physical pion mass 
predicted from the 5-parameter fit. 
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M n [MeV] M„ [MeV] 

FIG. 4: Predicted pion mass dependence of the Dir{I = 1/2) and DK(I — 0) scattering lengths using 
parameters from the 5-parameter fit. The solid curves are calculated using the parameters from the best fit, 
and the bands reflect the uncertainties. 

coupled-channel D s rj, and partly due to the energy dependence in the interaction. Using Eq. (19), 
the value of Z is in the range [0.25,0.41]. Therefore, the main component of the -D* (2317) is the 
isoscalar DK molecule. 

The result for the Dtt(I = 1/2) channel at the physical pion mass, with much smaller un- 
certainty, is still consistent with the indirect extraction, (0.41 ± 0.06) fm, in Ref. [61]. However, 
our prediction at a pion mass of about 266 MeV, (2.17^25) * s l ar g er than (0.81 ± 0.14) fm 
obtained in [14]. From Fig. 4, one sees that such a pion mass is close to the transition point where 
the scattering length changes sign due to the generation of a pole (for more discussions, see [49]). 
In such a region, the value of the scattering length changes quickly. For instance, decreasing the 
pion mass to 220 MeV, one would get a much smaller value (I.OO^q'qj) fm. 

From Table VII, one sees that all the dimensionless LECs are of natural size. The value of ho 
is much smaller than hi = 0.42. This is consistent with the N c counting because the ho term is 
suppressed by 1/N C as compared to the hi term. Furthermore, we also have the hierarchies /14 < h' 5 
and /12 < h% (recall that /12 = /124 — h' A and /13 = /135 — 2/ig). For both cases, the left hand sides 
are suppressed by 1/N C as compared to the right hand sides. This means the hadronic molecular 
assumption of the D* (2317) is consistent with the naturalness and iV c counting. 

All the above calculations have assumed the same mass for the up and down quarks, and ne- 
glected the electromagnetic interaction. This is the isospin symmetric case. However, the -D* (2317) 
decays into the isovector final state D s n. In order to calculate this isospin breaking decay width, 
one has to take into account both the up and down quark mass difference and virtual photons. This 
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has been done in Ref. [9]. In Ref. [9], the iV c -suppressed operators, i.e. the ho, hi and h± terms 
have been dropped, and a somewhat arbitrarily chosen natural range [—1,1] was taken for h' 5 . The 
isospin breaking decay width was calculated to be r(L>* (2317) -)• D s ir) = (180 ± 110) keV [9]. 
With the values of all the h^s in Table VII, the result is updated to be 

r(L>* (2317) -> D s tt) = (89 ± 27) keV. (20) 

We have used the isospin breaking quark mass ratio (m<2 — m u )/(m s — rh) = 0.0299 ± 0.018, where 
rh = (m u + TOd)/2, which is calculated using the lattice averages (up to end of 2011) of the light 
quark masses [69, 70]. 

VII. SUMMARY AND DISCUSSION 

The low-energy interaction between a light pseudoscalar meson and a charmed pseudoscalar 
meson was studied. We have calculated scattering lengths of five channels DK(I = 0), DK(I = 1), 
D S K, Dir(I = 3/2) and D s ir with four ensembles. Among these channels, the interaction of 
DK(I = 0) is attractive, and that of the others is repulsive. The interaction of D s ir is very 
weak, which is expected. The D s ir and DK(I = 1) channels are mixed since they have the same 
quantum numbers. To perform a more reliable analysis of these two channels, we need to construct 
the correlation matrix and use the variational method to extract the energies of the two channels. 
The chiral extrapolation was performed using SU(3) unitarized chiral perturbation theory, and the 
LECs h^s in the chiral Lagrangian were determined from a fit to the lattice results. With the 
same set of parameters, we made predictions on other channels including DK(I = 0), DK(I = 1), 
Dtt(I = 1/2) and D S K. In particular, we found that the attractive interaction in the DK{I = 0) 
channel is strong enough so that a pole is generated in the unitarized scattering amplitude. The 
pole can be a bound state pole. By further fixing the pole in the DK(I = 0) channel to the 
D* (2317), we found that the main component of the D* (2317) is the isoscalar DK molecule. We 
revisited the isospin breaking decay width of the L>* (2317) ->■ D s n. The result (89 ± 27) keV 
updates the old result (180 ± 110) keV obtained in Ref. [9]. It is nice to see that the uncertainty 
of the width shrinks a lot. 

It is possible to further constrain the values of h^s once simulations in other channels are done. 
Although a precise calculation of the other channels requires disconnected diagrams, one may 
obtain valuable information from the connected part only. The connected and disconnected parts 
can be calculated separately using partially quenched chiral perturbation theory (for reviews, see 
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Refs. [71, 72]), then a fit to the lattice calculation can be performed. This point has already been 
stressed, for instance, in Ref. [73] for the hadronic vacuum polarization and in Ref. [74] for the 
scalar form factor of the pion. 

In our chiral extrapolation, the resummed chiral amplitude is of 0(p 2 ). At this order, there 
is no counterterm to absorb the divergence of the loop function G(s), because loops only start 
from 0(p 3 ). As a result, we had to regularize the divergent loop by a subtraction constant, the 
pion-mass dependence of which was neglected. Were a full one-loop calculation available, the 
chiral amplitudes can be renormalized at one-loop order, and the representation of the pion mass 
dependence would be improved. However, more unknown LECs will be introduced in this way, and 
it is difficult to perform a fit with all of them to the present data. As mentioned above, more data 
can come from calculating the other channels, which is useful even if the disconnected contribution 
is neglected. Such a study with an improved chiral extrapolation is relegated to the future. 
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